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Abstract. The optical Faraday effect describes the rotation of linear polarization 
upon propagation through a medium in the presence of a longitudinal magnetic field. 
The effect arises from a different phase delay between the right and left handed 
polarisation components of the light. In this paper we predict a Faraday effect for 
a completely different system: electron vortices. Free electron vortex states were 
recently observed in transmission electron microscopy experiments, and they introduce 
new degrees of freedom into the probing of matter with electron beams. We associate 
a rotation of a vortex superposition with the fact that different phases are acquired by 
oppositely handed vortices propagating in a magnetic field. We show that, in contrast 
to the optical Faraday effect, the rotation of the electron beam occurs in vacuum and 
arises from the intrinsic chirality of the constituent vortex states. 
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1. The optical Faraday effect and its generalisation for electron waves 

Michael Faraday reported in 1845 that the polarisation of light can be affected by 
magnetic fields, an effect that now bears his name. Since then, the Faraday effect has 
found numerous metrological and research applications, including the ultra-sensitive 
detection of magnetic fields, [HE], or of fields generated by electron plasmas in 
interstellar space and the ionosphere [Sj H] . 

Faraday noted that the polarisation direction of light is rotated after passing 
through "heavy glass" exposed to a longitudinal magnetic field. We now understand 
that the Faraday effect arises from the different speed of propagation of right and left 
handed circularly polarised light through an optically active medium. The associated 
difference in accumulated phase between the circular components of linearly polarised 
light results in a rotation of the polarisation direction. 




Figure 1. (a) Illustration of the Faraday effect for spin angular momentum, i.e. optical 
polarization, where the two opposite circular polarization components propagate at 
different speeds through an optically active medium in a magnetic field. As a result, 
the linear polarisation rotates by an angle proportional to the magnetic field, (b) 
The analogous effect for transverse beam profiles of electrons, where the different 
propagation speed of states with opposite 0AM in a longitudinal magnetic field leads 
to image rotation. 
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One of the intriguing properties of light is that it can carry angular momentum: 
a spin contribution associated with circular polarisation (±/i), but also orbital angular 
momentum (0AM) [HIE]. While circular polarisation describes a rotation of the electric 
field vector upon propagation, the orbital angular momentum (0AM) is a feature of 
'twisted' light beams. The 0AM can take on arbitrary multiples of h depending on how 
tightly wound the phase fronts are. These so-called "vortex beams" have a rotational 
intensity pattern and are associated with a phase dependence exp{il(f)), where I is a 
non-zero integer and the azimuthal angle. 

Strictly speaking, Faraday rotation is not a relevant concept for optical 0AM. The 
reason is that there is no intrinsic mechanism in a gyromagnetic medium to produce 
the required 0AM state dependent dispersion, because selection rules forbid coupling of 
the 0AM to the atomic electron degrees of freedom. This is consistent with results from 
a recent experiment in which no rotation was observed for a superposition of right and 
left handed 0AM states (a Hermite-Gauss mode) propagating through cholesteric liquid 
crystals [7]. We note that a relative phase shift between right and left handed 0AM 
components will appear as a rotation of the intensity pattern [8J. Such phase shifts 
can be induced by spinning the medium through which the light propagates, inducing 
a 'mechanical' Faraday rotation, as demonstrated recently in a slow light medium [9]. 

Electron vortices are unusual quantum states that have only recently been 
predicted and produced in transmission electron microscopy (TEM) experiments 
pT| [T2] . Electron vortex beams have the same geometrical properties as their optical 
counterparts, being characterised by an exp(i/0) angular dependence related to Ih units 
of 0AM, but they also produce features that have no analogue in optics. In particular 
the circulation of charge in an electron vortex beam gives rise to an arbitrarily large 
orbital magnetic moment (Fig. |2]), distinct from the magnetic moment due to spin 
[TOj [T3] . Hence electron vortices can couple to electronic degrees of freedom through 
dipole selection rules forbidden to optical vortices [H]. 

Given the analogies (and differences) between optical and electron vortices, the 
question arises: Do electron vortex waves undergo something analogous to an optical 
Faraday effect? Here we show that there is indeed a Faraday rotation arising through 
Zeeman interaction from propagation parallel to a uniform, external magnetic field (i.e. 
in a geometry where there is no Lorentz force). 

2. Electron vortex states 

The dynamics of a non-relativistic electron propagating in a uniform longitudinal 
magnetic field B (with associated vector potential A) is described by the Hamiltonian 



where m and e = — |e| are the electron mass and charge, respectively. The Hamiltonian 
contains the kinetic energy due to the canonical momentum and the Zeeman interaction 




(1) 
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of the electron spin with an arbitrary external magnetic field B. The electron's two- 
component spinor wavefunction -ip satisfies the Pauli equation Hip = ih-^ip. Here 
fis = —QI^B^/'^ is the operator for the magnetic moment, hb = h\e\/{2m) is the Bohr 
magneton, g ^ 2 is the Lande g-factor for electron spin and & is the vector of Pauli 
spin matrices {dx, d'y, o"^). In agreement with typical parameters in transmission electron 
microscopes (~ 8 nA/nm^), we assume that beam currents are sufficiently low so that 
Coulomb repulsion can be neglected. 

In the case of the uniform magnetic field directed along the suitable choice 

for the vector potential is A = {Bzr/2)(j). Exploiting the cylindrical symmetry of the 
system, the Hamiltonian can be put in the form 

Here ul = \e\Bz/{2m) denotes the Larmor frequency, Lz = — and Sz = shcTz are 
the operators for the z component of 0AM and spin, respectively, where s is the spin 
quantum number. The first term gives the kinetic energy of motion along z, which is the 
same as in field-free space; the second and third terms together give the energy for the 
transverse motion, and have the form of the Hamiltonian for a harmonic oscillator with 
characteristic frequency u^, the final term gives the Zeeman energy, with contributions 
from both 0AM and spin. 




Figure 2. The circulating charge in an electron vortex beam generates a magnetic 
moment which interacts with external magnetic fields. The electron density function 
vanishes at the vortex core and the propagation direction is twisted around this core. 
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With Hamiltonian ([2]) the spinor components decouple and we can find 
monochromatic wave solutions which obey the time independent Schrodinger equation 
for a scalar wavefunction ip. This problem can be solved exactly for eigenstates with 
given z components of 0AM and momentum [15]. Separating the degrees of freedom, 

tpnisir, 0, z) = Rn\i\{r) exp(z/0) exp{iknisz) (3) 

we identify the radial modes 

= ,r^Z± ( ^] "' e--'<US (K) . (4) 

"'''^ ^ ]j TT{n+\l\)\WB \wb J WbJ 

where wb = '2'\fhJ\eB]\ is a characteristic width which depends on the magnitude 
of the magnetic field, and Ln is an associated Laguerre polynominal. The radial 
modes are characterised by the 0AM quantum number, /, and the radial mode number 
n = 0, 1, 2, . . . which denotes the number of radial nodes of the electron density function. 
The radial profile of a mode with / = 1 and n = is shown in Fig. 2. The transverse 
beam profile, given by Rn\i\{r^ exp(z/0), is the same as that of the Laguerre- Gauss beams 
familar from optical vortices, as was also pointed out in [T6] . 

Recognising that the combined second and third term in ([2]) describe a 2D harmonic 
oscillator, which in polar coordinates has energy eigenvalues h\uji\{fln + |/| + 1), the 
eigenvalues for the total energy are the well known Landau levels 

E = + h\uL\i2n + \l\ + 1) + huLil + gs). (5) 

The corresponding allowed wave numbers are then 

knis = ko^ 1 - ^[{2n + \l\ + l)h\ujL\ + {I + gs)hujL], (6) 

where E is the total energy determined by the electron source, and ko = \/2mE/h is 
the wave number for a plane (non- vortex) wave propagating freely along the z axis. For 
electron vortices, the phase aquired upon propagation depends on the magnetic field. 
We see from ([6]) that there is a path that depends upon the direction of the angular 
momentum through the signs of / and s. This corresponds to a phase Onis = knisZ 
acquired upon propagation which depends (via the Larmor frequency) on the magnetic 
field and is proportional to the total angular momentum comprising both a spin and an 
orbital component. 

If the magnetic energy is small compared to the total energy (a situation justified 
e.g. for electrons in TEM experiments), i.e. {2n + |/| + l)^|wi| + (/ + gs)hjjL <^ -E", 
we can apply the paraxial approximation to the wave numbers ([6]). The corresponding 
phase shift accumulated along the trajectory of the vortex then comprises three parts: 

Onis = koz - i2n + |/| + l)\ki^\z - (/ + gs)ki^z, (7) 

where ki = f^f^sBz/h is a spatial frequency which corresponds to the Larmor temporal 
frequency u^. The first term describes the phase evolution in free space; the second 
depends on the energy of the transverse motion due to the magnetic field, and the third 
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term arises from the Zeeman interaction with the total angular momentum. In the 
following it is this latter term which is important, as it causes a different phase shift for 
vortex states with opposite helicity. We note that electrons in superpositions of left and 
right spin components should result in a spin Faraday rotation, analgous to optics. Here 
we concentrate on 0AM Faraday rotation for which there is no optical counterpart. 

Just like for the optical Faraday effect, the differential phase shifts become 
observable as a rotation angle for electrons in superpositions of vortex states with 
opposite handedness. As they originate from the interaction of the magnetic dipole 
moment with the external field, the electron Faraday effect does not require the 
mediation of an optically active medium but occurs in vacuum! 



3. Considerations on observing the Faraday effect for electrons 

Electron vortex states have recently been generated in transmission electron microscopes 
(TEM) via diffraction from nano-fabricated holograms [121 113 HI H]- Using 
suitably designed holograms, also superpositions of vortex states can be generated (see 



Appendix B ). The required shape of the holograms is determined by the interference 
pattern of the target state with a reference wave function, e.g. a plane or spherical wave, 
resulting in transverse or longitudinal separation of the diffraction orders respectively. 

In order to realise electron Faraday rotation we require electrons in a superposition 
of two modes with the same spin and radial mode number but opposite vorticity ±1. 
The probability density then has an azimuthal dependence 

1 |2 

- I^n/s + ^n(-/)s| OC COS^[/(0 - <I>b)], (8) 

where we define 



1 / TTl 

<l>B = hz = ^J^f^BB,z. (9) 

For n = this is a petal pattern consisting of 2|/| maxima equally spaced around 
a circle, which after propagating through a region of a longitudinal magnetic field 
is rotated through the angle $b- The maxima are separated by phase singularity 
lines, where the phase changes by vr and the probability density vanishes. For the 
case of / = ±1 the transverse profile, shown in Fig. [3l is that of the HGio Hermite- 
Gaussian mode. Here the analogy with optical polarization is clearest, with the / = ±1 
components corresponding to the right and left handed circular polarization states, and 
the nodal line to the linear polarization. 

While the phase change depends on /, the rotation of the intensity pattern is 
independent of /. The Laguerre-Gauss modes form a complete basis with which an 
arbitrary wavefunction can be described, and therefore any intensity and phase profile 
will rotate through the same angle $b [H]- 

So far we have considered the eigenstates of ([2]) with a transverse scale determined 
by the magnetic field strength through the parameter wb- A beam with the same radial 
profile but a width w ^ wb is, however, no longer an eigenstate and will therefore 
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Figure 3. The electron density distribution rotates when propagating through a 
parallel magnetic field, and at the same time varies its beam waist periodically, with 
twice the period of the rotation (see inset). 



change upon propagation. Solving the paraxial wave equation, (see [Appendix A] ) we 
find that the radial profile retains the same Laguerre- Gauss form, only now expands 
and contracts periodically in time, or equivalently with propagation distance ( IA.7I) . 
This contraction happens at twice the Larmor frequency and hence twice every full 
rotation. The width variation can be understood in terms of the competition between 
diffraction, which dominates when w < wb, and the confining effect of the magnetic 
field, which dominates when w > wb- This is illustrated in Figure 3. 

The rotation angle $b, depends on the initial kinetic energy E of the electrons, as 
detailed in (|9]); as slow electrons spend more time in the magnetic field, their rotation 
angle is larger. While optical Faraday rotation is characterised by the Verdet constant 
(a proportionality constant of rotation angle per propagation distance and magnetic 
field strength), the electron 'Verdet' parameter varies with kinetic energy (see Fig. HI). 

Even with low-energy transmission electron microscopy, measuring the proposed 
deflection due to an interaction with a perpendicularly-magnetised sample remains 
challenging, in particular the neccessity to distinguish Faraday rotation from the usual 
cycloid motion of an electron beam within magnetic lenses. Ignoring relativistic 
corrections, typical values {E = 60keV, sample thickness = lOOnm, longitudinal field 
B = IT) yields a rotation about the beam axis of 0.06mrad. In certain geometries 
differential phase contrast techniques routinely measure deflections of this magnitude, 
when the rotation can be projected by long camera lengths to give measurable 
deflections. A more promising experimental approach will be to consider low energy 
electron beams. Low energy photoelectrons ejected by circularly polarised photons are 
also known to carry orbital angular momentum [20] and the subsequent propagation of 
these electrons through magnetic thin fllms may also visualise Faraday effects. 



Vacuum Faraday effect for electrons 
4. Conclusions 



8 



We have demonstrated that for electron vortex states propagating in a longitudinal 
magnetic field the Zeeman interaction produces an 0AM dependent dispersion. This 
results in a rotation of the probability density of a superposition of vortex states about 
the beam axis, analogous to the optical Faraday effect. To the best of our knowledge 
this is a new concept, it is an effect that is not present in optics, and it may lead 
to applications in electron spectroscopy. The magnitude of the rotation scales with 
the magnitude of the magnetic field, and increases with decreasing energy. There 
are a number of interesting applications which may follow. 0AM Faraday rotation 
provides the possibility of spatially resolved measurements of longitudinal magnetic field 
components, analogous to the measurement of transverse fields in Lorentz microscopy, 
by measuring the rotation angle of a vortex superposition. Moreover, we note that in 
the approximation considered here, spin and 0AM are separately conserved. This would 
not be expected for relativistic non-paraxial beams [I3] or spatially varying magnetic 
fields [21] , suggesting a route to investigating intrinsic spin-orbit coupling in an electron 
vortex beam, via a spin-dependence of the rotation angle. 
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Figure 4. Verdet parameter as a function of energy, i.e. the rotation angle per T and 
m. 
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Appendix A. Calculation of beam width variation using the paraxial 
approximation 

To find the propagation of a beam with a given beam waist we evaluate the time- 
independent Schrodinger equation with Hamiltonian ([2]) as before in the paraxial 
approximation. We assume a solution of the form 

^{r,(f),z) = u{r,(j),z)e'''°\ (A.l) 

where u{r, (p, z) is an envolope function which describes the evolution of the beam profile 
upon propagation, ko is as defined in the main text. If u varies sufficiently slowly with 
z we can use the paraxial approximation 









dz^ 


< 


oz 



\V'u\ 



(A2) 



Then, substituting f lA.ip into the Schrodinger equation, we arrive at the paraxial wave 
equation 



u 



2ikcS-— 



k'lk^r'^u 



2kokL 



L,u = 0. 



'dz h ^^-^^ 

The first two terms here are the same as in the paraxial equation for an optical beam 
propagating in vacuum, or for an electron beam in field-free space. The third term 
represents the confining effect of the magnetic ffeld, and the final term gives the Zeeman 
interaction. For an 0AM eigenstate with u oc e*'*^, the last term takes a constant value 
— We can then factor out the phase due to the Zeeman interaction by writing 

u = ve-'^^'-^ (A.4) 

for some function v{r, (j),z). v then satisfies the equation 



dv 

Vlv + 2iko— - klklr'^v 
oz 



0. 



(A.5) 



This equation has the same form as the equation from paraxial optics for a medium with 
quadratically varying refractive index n{r) = hq — ^n2r^, where Uq and n2 are constants 
(see for example [22]), only here ^2 — )■ Such an optical system supports Laguerre- 
Gauss type modes which experience a periodic width variation due to the competition 
between diffraction and the focusing effect produced by the refractive index variation. 
A similar effect is described in [23], for the propagation of vortex wavepackets in a 
transverse field. 

The solutions can be written 1241 [25] as 



Vni{r,(f),z) 




(A.6) 




Figure Bl. (a) Hologram grating to produce a superposition of vortex modes with 
I — ±1. (b) Grating to produce the same superposition but with the diffraction orders 
separated longitudinally along the beam axis, (c) Diffraction pattern produced by 
grating in (a), with the desired superpositon in the ±1 diffraction orders. 



where w{z) is the beam width, R{z) is the wavefront radius of curvature and $,{z) gives 
the longitudinal phase shift. The equation (1A.6I1 is the same as for the LG modes in free 
space, except here the functions w{z), R{z) and C,{z) are different. Choosing z = to 
coinside with one of the minima of w{z), and calling this minimum value Wq, the width 
function can be written, for the electron beam in a magnetic field, as 



w{z) = Wb 



Wq 
Wb 



COS {2kLz). 



(A.7) 



Appendix B. Diffraction grating patterns for production of electron vortex 
superpositions 

For interference with a plane wave ip oc e*'^^^, in which case the diffraction orders will 
be separated transversely, a suitable hologram pattern is generated from 

1^2, ^fl if||2cos/(0-0o)+e^'^-p>i 
°° I otherwise, 

where 0o specifies the orientation of the singularity lines. A grating producing a 
superposition of / = 1 and / = — 1 modes in the first diffraction order is shown in 



Fig. 1(a) , where we have chosen 0o = 0. Note that the left and right sides of the grating 
are displaced by half a period with respect to one another, introducing the necessary vr 
phase shift. In a similar way a hologram which separates the diffracted beam components 
longitudinally rather than transversely could be used [HI [19] . An example of this type 
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of grating is shown in Fig. 1(b) This is calculated using the same method, only with 
the plane wave factor ip oc e*'^^^ replaced with the spherical wave profile ip oc e**^*^^, 
where C is a constant which determines the curvature of the wavefront, and hence the 
spacing of the diffraction orders. Such grating patterns can be produced in the same 
way as those already used to generate single vortex modes. We note that gratings that 
generate superpositions of vortex beams with beams that have flat phasefronts have 
recently been generated and used to investigate the Gouy phase for electrons [26] . 
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